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SUMMARY 


The development of a potential- flow/boundary-layer method for calculating 
subsonic and transonic turbulent flow past airfoils with trailing-edge separa- 
tion is reported. A mcment-of-mcmentum integral boundary-layer method is used 
which employs the law-of- the- wall/law-of- the- wake velocity profile and a two- 
layer eddy- viscosity model and ignores the laminar sublayer. All integrals 
across the boundary layer are obtained in closed form. Separation is assumed 
to occur when the shearing-stress velocity vanishes. A closed-form solution 
is derived for separated- flow regions where the shearing stress is negligible. 

In the potential- flow method, the exact form of the airfoil boundary condition 
is used, but it is applied at the chord line rather than at the airfoil surface. 
This allows the accurate computation of flow about airfoils at large angles of 
attack but permits the use of body-oriented Cartesian computational grids. The 
governing equation for the perturbation velocity potential contains several 
terms in addition to the classical small-disturbance terms. 


INTRODUCTION 

This paper reports on the development of a potential-flow/boundary-layer 
method for calculating two-dimensional subsonic and transonic separated flow 
past airfoils. The form of separation which is discussed is the turbulent 
trailing-edge type. Previous potential- flow/boundary-layer treatments of 
the trailing-edge separation problem (refs. 1 to 6, for example) have been 
restricted to low-speed flow. Previous treatments of the transonic trailing- 
edge separation problem (ref. 7, for example) have used solutions to the Navier- 
Stokes equations, which are very expensive to compute. 

The present method was developed to show that it is possible and feasible 
to calculate separated flow past relatively thick, cambered airfoils with angles 
of attack near stall and fairly large subsonic free-stream Mach numbers up to 
the drag-divergence Mach number. These are the conditions under which many 
rotorcraft and aircraft airfoils operate. It is not the intent at this time 
to discuss leading-edge separation or shock-induced separation (separation at 
the shock wave). The basic conclusion of the present report is that it is, in 
fact, possible and feasible to calculate transonic flow past airfoils having 
trailing-edge separation with a potential- flow/ boundary- layer method. 

The boundary-layer and the potential- flow methods which are used are 
described in the next two sections. The methods are coupled in that the 
boundary-layer method uses edge- velocity distributions computed from the 
potential- flow program, and the potential- flow program uses a displacement 
thickness and a separation-point location determined by the boundary-layer 
method. 



SYMBOLS 


A,B,C coefficients of potential- flow governing equation 

a speed of sound 

a ij,bi, c i coefficients of boundary-layer equations (i = 1,2,3 and 
j = 1,2,3) 

B law-of- the- wall boundary-layer coefficient 

Cj lift coefficient 

Cp pressure coefficient 

Cp* sonic pressure coefficient 

c airfoil chord 

C3*,C3 normalized forms of C3 for u* = 0 

d^ coefficients in equation for gradient of displacement thickness 

(i = 1,2, 3, 4) 

I I5 definite integrals 

K coefficient in Clauser model for eddy viscosity 

M Mach number 

p pressure 

Rq transition point Reynolds number based on momentum thickness 

S Sutherland temperature parameter 

U,W transformed velocity components 

u,w x and z velocity components 

u*,u$ law-of-the-wall and law-of- the- wake velocity parameters 

v perturbation velocity parameter 

v ro free-stream speed 

x,z Cartesian coordinates for chordwise and normal directions in 

inviscid treatment and for tangential and normal directions in 
boundary-layer treatment 

a angle of attack 
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edge-velocity gradient parameter 
circulation 

ratio of specific heats 

intermittency factor 

increment added to a 33 

increment added to c 3 

mesh spacings in x and z directions 

boundary-layer thickness in transform plane (distance from 
surface at which U = U e ) 

boundary- layer thickness parameter in physical plane 
boundary-layer displacement thickness 
eddy viscosity 

eddy viscosity in inner layer 
eddy viscosity in outer layer 

normal and tangential boundary-layer coordinates in transform plane 

value of n where U = 0 

boundary-layer momentum thickness 

law-of- the- wall boundary-layer coefficient 

turbulent viscosity 

kinematic viscosity 

coefficient used in starting solution to relate u* and t£ 
dens ity 

density at wall 

single-variable velocity potential in transformed plane 
single-variable velocity potential in physical plane 
perturbation-velocity potential 
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Subscripts: 

e boundary- layer edge 

i interface between inner and outer layers of boundary layer 

s separation point 

te trailing edge 

co free stream 

± upper and lower surfaces of airfoil 


INTEGRAL BOUNDARY-LAYER METHOD FOR TURBULENT FLOW 
Derivation of Equations 

The present boundary-layer method is a simplified form of the turbulent- 
flow integral method of Kuhn and Nielsen (ref. 8) , which is a refinement of the 
method of Nash and Hicks (ref. 9) . The equations which are solved are the 
momentum equation, the moment- of-manentum equation, and a boundary- layer edge 
condition. The present model, like the model of reference 9, does not include 
the laminar sublayer. The wall is assumed to be adiabatic and the turbulent 
Prandtl number to be 1 so that the total enthalpy is constant across the boundary 
layer. The Reynolds stress is approximated with an eddy- viscosity model which 
is like the Cebeci-Smith model (ref. 10) but does not contain laminar-sublayer 
terms. Compressibility is accounted for with the Stewartson transformation 
(ref. 11) as in reference 8. 

The fact that the laminar sublayer can be ignored in the present integral 
treatment of separated flow may appear unusual to researchers acquainted with 
finite-difference solutions to the boundary- layer and Navier-Stokes equations. 

The influence of the sublayer on the physical solution diminishes as separation 
is approached because the sublayer thickness decreases whereas the boundary- 
layer, displacement, and momentum thicknesses increase rapidly. As a result, 
the physical error incurred by neglecting the laminar sublayer at separation 
is actually less than that incurred in the attached- flow region. However, for 
finite- difference treatments of the separation problem, the influence of the 
sublayer on the numerical solution increases drastically because seme minimum 
number of grid points must be maintained within the shrinking sublayer at the 
same time that the boundary layer as a whole is increasing rapidly in thick- 
ness. Thus, in finite-difference solutions of the separation problem, resolu- 
tion must be increased in a region of diminishing physical importance at the 
same time that higher resolution is needed for physical reasons in other parts 
of the solution domain. 

The continuity, momentum, and energy equations which govern the turbulent 
compressible flow under consideration are 
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where x and z are the tangential and normal coordinates, u and w are 
the respective velocity components, p, p, and y are the pressure, the 
density, and the ratio of specific heats, and is the turbulent viscosity. 

The pressure, density, and speed in the free stream are designated as p , p , 
and v . “ 00 

OO 

The transformed coordinates for the present version of the Stewartson 
transformation are 



(4) 
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where a is the speed of sound and the subscript e designates boundary-layer 
edge conditions. The transformed velocity components U and W are 
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and 
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The governing partial differential equations are written in transform variables 
as 
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is the eddy viscosity. 

Coles' law-of-the-wall/law-of-the-wake velocity profile (ref. 12), which is 
used in the present method, can be written as 


U 



+ ug sin 2 
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where 6 is the boundary-layer thickness (distance from the wall at which 
U = U e ) , u* and ug are the wall and wake velocity parameters, Voo is the 
free-stream kinematic viscosity, and k and B are constants which are 
given the values 0.41 and 5.00, respectively, in this report. Note that at the 
boundary-layer edge, equation (11) has the form 



+ u 3 
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and that the derivative of equation (12) with respect to £ can be written as 
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where 
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Note also that the boundary-layer displacement thickness 6* and momentum 
thickness 0 can be written as 


s * . r z=s °/'i . — 

^■r = n \ P^U Q 


P U \ Poo a co j U* U 3 

dz = S{ + 1 + (Y - 1)M 2 

>e u e/ Pe a e ] < u e 2u e 


Y - 1 ^2u* 2 3u 3 2 0 + 1 1 ) (u*up) 

Mg 2 + + 

2 K 2 U e 2 8U e 2 KU e 2 


pz=6 c PU / u\ Poo a oo ju* u 3 

9 = \ [1 - — ] dz = 6 / + 

z=0 Pe u e V u e/ Pe a e j^e 2 ^e 


2u* 2 3 u 3 2 (1 + Ii) (u*u 3 ) 

K 2 U e 2 8U e 2 KU e 2 


where S c is the transformed value of 6 in the physical plane, M is the 
Mach number, and the definite integral 1-j is 
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The factor P 00 a 00 /Pe a e can be evaluated with the equation 
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A two-layer eddy-viscosity model is used which is similar to the Cebeci- 
Smith model (ref. 10) but which does not have viscous sublayer terms. In the 
inner layer a Prandtl mixing-length model of the form 
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is used. In the outer layer the Clauser model, which can be written as 
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is employed, where K is a constant which, in the present method, is given the 
usual value of 0.016, and y is the intermittency factor. 

In this report the intermittency factor is approximated with two line 
segments as 
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This simple approximation is compared with the experimental data of Klebanoff 
( re f* 13) and with the more elaborate empirical expressions of Klebanoff and 
of Cebeci-Smith (ref. 10) in figure 1. 


As in most two— layer models , the switch from the inner to the outer— layer 
form is made at the point where the values of two forms are equal, that is, 
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In general, this point is located in the interior of the boundary layer where 
the intermittency factor is 1 and the density is approximately equal to the edge 
value. With equations (11), (21), and (22), equation (24) can be written as 
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where is the outer limit of the inner layer. 

In the present method the differential equations which are solved are the 
integral forms of the momentum equation and the moment-of-momentum equation. 
With equations (8) and (11) the integral form of equation (9), the momentum 
equation, can be written as 
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where P w is the density at the wall. The ratio P w /Poo is given by the 
equation 
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and where the definite integral 13 is 
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The treatment of the shearing-stress integral in the term c 3 is one of 
the distinguishing features between the present method and the preceding methods 
of Kuhn and Nielsen (ref. 8) and Nash and Hicks (ref. 9) . In reference 9 the 
term was approximated as a perturbation about the value for equilibrium flow. 

In reference 8 the integrals across the boundary layer were evaluated numerically 
because of the complexity caused by the inclusion of the laminar sublayer terms. 
In the present method the shearing-stress integral is integrated in closed form. 

An additional approximation is made in order to evaluate the integrals in 
equation (38). The ratio P/P m is 
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In the present method it is assumed that for subsonic and transonic Mach numbers 
the ratios P/P m and P^jpJ^ can be approximated as 



(41) 
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With equations (11), (23), (41), and (42), the term c 3 in equation (38) is 





(43) 


where U* is the value of U at the point n = hi- In regions where the flow 
is separated the term 
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must be added to equation (43) to account for the absolute value in the first 
term on the right side of equation (38) . The point at which U vanishes is 
designated ri j - It has been assumed that rij is less than Hi* 

The quantities £ and U e are more difficult to deal with than x and 
u e . Therefore, equations (13), (26), and (33) are transformed into 
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for purposes of numerical integration. The quantities U e , M e , dU e /du e , and 
a e/ a oo are 9i- ven by the following equations: 
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If | a i j I , the determinant of the coefficient matrix for equations (45), 
vanishes, the derivatives du*/dx, di^/dx, and dS/dx become unbounded. This 
type of behavior is expected for a direct boundary-layer solution at a separa- 
tion point, but it is not expected before separation occurs. However, because 
of the approximate nature of the present treatment, the determinant for the 
present method vanishes at a point just upstream of separation as well as at 
the separation point. In order to improve the behavior of the numerical solu- 
tion, the moment-of-momentum equation is manipulated so that the resulting coef- 
ficient matrix does not vanish until the separation point is reached. The 
manipulation consists of adding the quantities 
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and 
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to equations (36) and (43). The effect of this manipulation is to lag the 
solution for the derivative d6/dx. 


Starting Solution 

In order to start the numerical integration of the boundary-layer equations 
it is necessary to obtain the initial solution for u*, ug, and 6 at the 
transition point. This solution is determined by iterative solution of equa- 
tions (12), (18), and the equation 
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for a specified value of the transition-point Reynolds number based on momentum 
thickness Rq. An expression for the quantity II is obtained from reference 14 
as 
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where 
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The momentum thickness is related to the Reynolds number based on momentum 
thickness by the equation 


e = v e R 0 /u e 


(55) 


where 


v e 
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where and T e are the free-stream and boundary-layer-edge temperatures, 

is the free-stream kinematic viscosity, and S is the Southerland temper- 
ature parameter. The transition-point Reynolds number based on momentum thick- 
ness is generally given a value of 320. 


Integration of Equations in Attached-Flow Region 

The solution for the quantities u*, ug, and 6 is advanced downstream 
in the attached-flow region by numerical integration of equations (45) . A 
second-order predictor-corrector technique is used in the integration process. 
The edge velocity u e and the gradient of the edge velocity are obtained from 
the current inviscid solution. The numerical integration is terminated at the 
separation point, where u* vanishes. 

It it interesting to note that the differential equations do not depend 
explicitly on the kinematic viscosity. Therefore, the Reynolds number influ- 
ences only the starting solution explicitly. This result is a consequence of 
ignoring the laminar sublayer. 
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In the present method there are two options for evaluating the law-of-the- 
wake velocity parameter ug. The most obvious procedure is to use the value 
obtained from the integration of equations (45) . The second procedure is to 
use equation (12), which is the integrated form of the first of equations (45). 
The second procedure reduces errors incurred in numerical integration which tend 
to be disruptive as the separation point is approached and the determinant of 
the coefficient matrix ai j | approaches zero. Another advantage to the use 
of equation (12) is that itmakes the solution explicitly dependent on the kine- 
matic viscosity and hence the Reynolds number. 

In regions where the inviscid flow is expanding, the parameter ug 
approaches zero but should not become negative. If the numerical procedure 
discussed previously predicts a negative value for ug, the quantities ug and 
dug/dx in the second and third of equations (45) are given values of zero, and 
these equations are solved for du*/dx and dS/dx and hence for u and 6. 


Integration of Equations in Separated-Flow Region 

It has been pointed out in references 8 and 9 that numerical solutions to 
boundary-layer equations posed in the form of equations (45) are not stable in 
regions of separated flow. However, stable solutions can be obtained if the 
law-of- the- wall velocity parameter u* is prescribed and equations (45) are 
solved for the quantities ug, 6, and u e (or U e ) . 

The present method is intended for use in treating flow about airfoils with 
trailing-edge separation. It is well-known empirically that the reverse-flow 
speed is small relative to free-stream speed in most trailing— edge separated 
flows. Therefore, it is assumed in the present method that u is zero in 
separated- flow regions. It follows from the first of equations (45) that du /dx 
is also zero in these regions. It is also assumed in the present method that the 
entire region from the first separation point detected to the trailing edge is 
separated. 

When u* and du*/dx are zero, the quantities ug and dug/dx are equal 
to U e and dU e /dx, and the second and third of equations (45) can be written 
as 
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The quantity C3 depends only weakly on U e and is considered constant in the 
present treatment. 

There is a simultaneous closed-form solution to equations (57) and (58) 
which can be written as follows: 


6 = 6 S [1 + c 3 *(x - x s ) ] 
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where 5 S , U e ^ s , and x s are separation-point values and where 
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The quantities 6 S and x s are determined from the boundary-layer solution. 
However, the quantity U 6fS is determined from the inviscid solution in order 
to couple the viscous and inviscid solutions and to facilitate convergence. 

In order to derive an expression for U e , s , a single-variable velocity- 
potential function $ is introduced which is related to the edge velocity in 
the transformed plane U e by the following equation: 


d$ 



With equation (61), equation (63) can be integrated to determine the function $. 
The following equation for U e>s is obtained from this function evaluated at 
the trailing edge: 


5 c 3* ($te " $s) 

6 [l + c 3 *(x te - x s )]5/6 _ 1 


(64) 
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where $ s is the separation-point value and $ te and x te are trailing-edge 
values. The value given by equation (64) and the value of U e at the separa- 
tion point converge as the inviscid solution converges. 

In order to obtain the solution for the edge velocity in the physical plane 
u e , this quantity is associated with a single- variable velocity-potential func- 
tion 9 , and it is noted that U e (and hence u e ) does not vary a great deal 
in separated- flow regions. It follows that 


dp dp 9te “ 9s 

— = — U p s 

dx d$ “ $ s 


(65) 


or 


5 c 3 *(cpte - 9 S ) 

6 {[l + c 3 *(x te - x s )] 5 /6 - 1} [1 + c 3 *(x - x s ) ] 1/6 


The evaluation of the velocity potentials <Pte and ip s is discussed in the 
section of the paper entitled "Approximate Potential-Flow Method." Equa- 
tion (66) is used to evaluate the quantity u e in equation (97) for the 
Dirichlet boundary conditions for the inviscid solution in the separated-flow 
region. 

It has been found that equation (66) tends to overpr edict the variation of 
u e for incompressible and low-speed conditions for which it is well-known 
empirically that u e is almost constant. Therefore, the present method con- 
tains an option to evaluate u e in the limit as c 3 * vanishes. This value is 


Uo = 


9te " 9s 
x te ' x s 


(67) 


Evaluation of Displacement-Thickness Gradient 

In the attached-flow region the inviscid solution is influenced by the 
boundary-layer solution through the quantity d5*/dx, the gradient of the 
displacemnt thickness. It is common practice to determine 6 from the 
boundary-layer solution and then obtain d6*/dx by numerical differentiation. 
In general, this procedure produces eratic variations of dS*/clx which must 
be smoothed. 

In the present method the gradient dS*/dx is calculated directly frcm 
the gradients of u*, ug, 6, and u e . The expression for this gradient, 
which is the derivative of equation (17) with respect to x, is 


18 


( 68 ) 


where 


d6* du* dug d6 du e 

= d-j + do + do — + d 4 

dx dx dx dx dx 


a e Pe 6 

a «Pco KU e 


d-j = (1 + (y - l)M e 2 


2 u* 0 + I l ) u 3 

K U e 2U e 


a eP e 6 

2V 


d 2 = — <1 + (Y - 1)M e 2 


1 - 


(1 + l q )u* 3ug 


KU e 4U e 


6* 

do = — 


dU, 


d 4 = 


e l(Y + 1) M o^e 6* 

d “ : ' _T_ , * £L l" „ e 2 


a P P 


e^e 


a ooPc 


r 

l u* u 3 

M 

i 

£ 

K u e 2 U e 



(69) 


(70) 


(71) 


(72) 


It is not necessary to smooth the values obtained from equation (68) with respect 
to x. However, it is necessary to relax the values of d5*/dx with respect to 
the inviscid-solution iteration variable near the trailing edge. Equation (68) 
is used to evaluate the gradient d<5*/dx j. n equations (95) and (96) for the 
Neumann boundary conditions for the inviscid solution in the attached-flow 
region. 


Alternate Form of Third Equation 

In the present treatment, as in that of Kuhn and Nielsen (ref. 8) and of 
Nash and Hicks (ref. 9) , the third equation is the moment-of-momentum equation. 
An alternate approach involves the use of the integral mean kinetic-energy 
equation, which can be written as 
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The definite integrals I 3 , I 4 , and I 5 are given in the following equations: 


1 C° =7T o 

l 3 = 1 + — In 2 ( 0/1 T) cos a da = 1.16713 

2 * 0 
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I 4 = ~ + J 1 + 
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dO = 1.39592 


(80) 
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do 


= 7.23979 


(81) 
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In the regions where the flow is reversed, the following quantity is added to 
equation (78) s 
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In order to improve the behavior of the numerical solution near separation 
points, the quantities 
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and Ac 3 given in equation (51) are added to equations (76) and (78) 
quantity c 3 * in the separated-flow solution is 


The 


C3 * = S.\ 

3 H a ” " eV u * =0 


(84) 


No other changes to the formulation are required. 
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APPROXIMATE POTENTIAL-FLOW METHOD 

The transonic potential— flow method is the type presented by Murman and 
Cole in reference 15. The potential function is determined by solving finite- 
difference approximations to the governing equation at the grid points by 
successive line-overrelaxation techniques. In the present method a body- 
oriented, nonuniform Cartesian grid is used. A schematic representation of 
this computational grid is given in figure 2. The present method is unusual 
in that the exact nonlinear form of the airfoil boundary condition is used, but 
the boundary condition is applied on the airfoil chord line rather than the 
actual airfoil surface. This approximate nonlinear boundary condition has the 
advantage of the linear boundary condition of small-disturbance theory in that 
elaborate mappings of the airfoil are not required in order to apply the boundary 
condition at the surface, and yet it does not erroneously force the stagnation 
point to always be located at the airfoil leading edge. It should be noted that 
the stagnation point needs to be free to move back from the leading edge at the 
relatively large angles of attack at which trailing— edge separation occurs. The 
separation point is modeled in the potential-flow method as a free-streamline 
separation. Neumann boundary conditions are used on the attached-flow part of 
the airfoil and Dirchlet boundary conditions are used where the flow is sepa- 
h^hed. present, the Dirchlet boundary conditions are constructed using 
either an approximate closed-form, separation-region solution or the well-founded 
empirical assumption that the pressure in the separation region is constant for 
nearly incompressible flows. 

In the present potential- flow method various forms of the governing equa- 
tion ranging from the full potential equation to the small-disturbance equation 
can be treated. In general, the governing equation can be written in the form 


3^cJ> 

A + B 

3x^ 


32 (|, 
3x 3z 


3 2(j> 

+ C = 0 

3z2 


(85) 


where x and z are the coordinates in the chordwise and normal directions 
and <}> is the perturbation velocity potential function. Very simple finite- 
difference expressions are used in the present method. The derivatives 
3 2 <|>/3x 3z and 3 2 (j>/3z 2 are always approximated with central-difference 
expressions. The derivative 3 2 <j>/3 x 2 i s approximated with a central-difference 
expression if the coefficient A is positive and a backward— difference expres— 
s *on ^ * s negative. Shock waves can be treated either conservatively or 
nonconservatively. When the conservative approach is used, the derivative 
3 2 <t>/3x 2 is approximated with the shock operator developed by Murman (ref. 16) 
and simplified by Barnwell (ref. 17). 

It is well-known that the classical small-disturbance form of the govern- 
ing equations lacks accuracy and that the full-potential equation is difficult 
ho stabilize when the flow past the airfoil is critical and the procedure 
described in the previous paragraph is used. The difficulty with stabilization 
is due to the fact that the boundary conditions are not applied at the surface 
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and the fact that rotated differencing (ref. 18) is not used. However, a stable 
intermediate form has been found during the present investigation which yields 
the same results as the full-potential-equation form at nearly critical Mach 
numbers. This intermediate form is written with the following coefficients: 



( 86 ) 


B = -2M„ cos a + 

CO 


w 

(sin a + — 

\ v_ 


(87) 


c = 1 - M co 2 sin2 a (88) 


where M^, a, and y are the free-stream Mach number, angle of attack, and 
ratio of specific heats and where u and w are the velocity components in 
the x and z directions. For comparison, the coefficients A, B, and C 
for the classical small-disturbance equation are 


= 1 - M 2 + (Y + 1)M 2 - 


B = 0 


C = 1 


(89) 

(90) 

(91) 


For the full potential equation the coefficients are 
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(94) 
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Note that in the present method (eqs. (86) to (88)), the dependence of A on w 
and the dependence of C on u and w have been ignored, and A, B, and C 
are correct as infinity is approached, A comparison between results of the 
present approximate potential method and those of the full potential method 
(ref. 18) is given in figure 3(a). Also shown are results calculated with the 
approximate method but with B = 0 (fig. 3(b)). it is clear from the results 
that the mixed derivative is an important contributor and should not be ignored. 

The approximate Neumann boundary conditions for the attached-flow region of 

the airfoil surface for the present method and the full— potential— equation method 
are 


8(J) 
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- v sm a 

oo 
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(95) 


where z — z + (x) are the equations for the upper and lower airfoil surfaces. 
As indicated, these equations are exact in form but are evaluated along the 
airfoil chord line. For comparison, the boundary conditions for the classical 
small-disturbance equation are 


9<j) 

3z 



z=±0 


(96) 


It should be noted that the present method and the full-potential-equation 
method both exhibit small-disturbance-like leading-edge singularity problems 
when grid points are located too near the leading edge. These problems occur 
because the boundary conditions are not applied at the airfoil surface. 

In regions where the flow is separated, Dirchlet boundary conditions are 
used. One of the forms which is used in the present method is 


r* s=x 

4> = 4>s + j ( u e “ v oo cos a ) ds (97) 

^s=x s 

where 4> s and x s are separation-point values and u e is the boundary-layer 
edge velocity. The present boundary-layer method has a closed-form expression 
for u e . Therefore, the integration in equation (97) is performed analytically. 

The circulation on the airfoil is evaluated with a Kutta condition regard- 
less of whether the flow at the trailing edge is attached or separated. This 
condition is applied in the form 


u +,te = u^te 


(98) 
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which can be written as 


<j) + (x te + Ax) - 4> + (x te - Ax) <M*te + Ax) - <fr_ (*te ~ Ax) 
2 Ax 2 Ax 


(99) 


where the subscript te denotes trailing-edge values. Equation (99) can be 
rewritten as 


r = 4> + (x t e “ Ax) - 4>_ (x te - Ax) 


(TOO) 


where T is the circulation. For attached flow the values of the perturbation 
velocity potential function in equation (TOO) are obtained by linear extrap- 
olation, and for separated flow the values are obtained from equation (97) . 

The surface pressure coefficients are calculated from the following 
equation: 


C D = v + - M 2 v 2 + — M 4 v 3 
P 4 oo 24 


( 101 ) 


where 


v = - — (2 cos a + 

M 

This form is accurate at transonic speeds and is a determinate form for incom- 
pressible flew. 


u \ w / w 

— 1 - — 2 sin a + — 


( 102 ) 


In the present method, the values of the perturbation velocity potential 
at the airfoil surfaces which are used to calculate the surface pressure coef- 
ficients are obtained with the following equation: 


3<j> 

(J> + (x) = <(>(x,±Az) ± Az — (x,±0) (103) 

3 z 

where the derivatives 3<|)/3z are evaluated with equations (95). This form of 
extrapolation is also used for the full-potential-equation method. However, 
for the small-disturbance form the surface perturbation velocity potentials are 
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evaluated with linear extrapolation. (It is interesting to note that neither 
extrapolation procedure seems to work for the other form of the governing 
equation . ) 

The single-variable velocity-potential function 9 , which is defined in 
equation (65) is related to the surface perturbation velocity potential func- 
tion <(> by the equation 


(f = v x cos a + 4> 


(104) 


This equation is used to determine the separation-point value 9 S and the 
trailing-edge value <fte which appear in equations (66) and (67). At the 
separation point <j> is evaluated with equation (103) for both the present 
method and the full-potential-equation method and with linear extrapolation 
for the small-disturbance method. At the trailing edge <j) is evaluated with 
linear extrapolation for all methods. 

It should be noted that the incompressible limit of the present 
perturbation-velocity-potential method for separated flows was presented in 
reference 19. The advantages of the present method are its accuracy and the 
ease with which boundary conditions can be applied. The present method should 
be useful for accurately calculating flow past three-dimensional configura- 
tions, for which it is extremely difficult to enforce boundary conditions at 
wing surfaces. 


RESULTS 

Results computed with the present method for the lift curves of the 
NACA 0012 airfoil at free-stream Mach numbers of 0.15, 0.30, and 0.50 and 
Reynolds numbers of 2.0 x 10 6 and 6.0 x 10 6 are compared with experimental 
results in figure 4. Experimental results for free-stream Mach numbers of 
0.15 and 0.30 were obtained in the Langley Low-Turbulence Pressure Tunnel (LTPT) 
by Charles L. Ladson (unpublished) , and results for free-stream Mach numbers 
of 0.30 and 0.50 were obtained in the Langley 8-Foot Transonic Pressure Tunnel 
(8-Ft TPT) by Charles D. Harris (unpublished). 

The results for the variation of the lift coefficient Cj with angle of 
attack at the lowest Mach number are essentially the same as results for incom- 
pressible flow. In the experiment, the boundary layer was tripped at points on 
the upper and lower surfaces 0.05c downstream of the leading edge, in the com- 
putation, the lower-surface boundary layer was tripped at 0.20c to avoid the 
numerical difficulties caused by starting the boundary layer ahead of the stag- 
nation point, and the upper-surface boundary layer was again tripped at 0.05c. 

The first computations at a free-stream Mach number of 0.30 were made with 
the boundary layer tripped at the same points discussed previously. These 
results, which were in sharp contrast to the experimental results at this Mach 
number and which are not shown in the figure, were almost the same as the incom- 
pressible results. An examination of the pressure distributions revealed a 
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shock wave at 0.025c behind the leading edge on the upper surface for large 
angles of attack. Because this shock wave was located ahead of the upper- 
surface transition point, no shock-wave/boundary-layer interaction could occur. 
When the upper-surface transition point was moved to the leading edge, the 
results indicated by circular symbols, which are in excellent agreement with 
experimental results, were obtained. These results demonstrate that the decrease 
in maximum lift coefficient with increasing free-stream Mach number is due to 
a shock-wave/boundary-layer interaction. The computational results calculated 
with a linear form of the present method are represented by square symbols in 
figure 4. These results show that linear methods will not predict accurately 
the effect of Mach number on the maximum lift coefficient. 

The results for a free-stream Mach number of 0.50 were calculated with the 
upper-surface transition point 0.05c behind the leading edge. These computa- 
tions, which are in good agreement with experiment, were successful because the 
shock wave was located downstream of the transition point. 

The pressure distribution for an NACA 0012 airfoil traveling at transonic 
speed is shown in figure 5. These results show a shock wave on the upper sur- 
face near the leading edge and separation at the trailing edge. 

Sample calculations have been made with the alternate form of the present 
method, in which the mean kinetic-energy integral rather than the moment-of- 
momentum integral is used as the third equation. These calculations show there 
is no appreciable difference in the results of the two formulations. 


CONCLUDING REMARKS 

It has been demonstrated that potential-flow/boundary-layer methods can 
be used to calculate transonic flow past airfoils which are experiencing 
trailing-edge flow separation. Such methods can be used to determine the 
maximum lift coefficient and perhaps the effect of flap deflection on airfoil 
characteristics. Also, a transonic potential-flow method with easy-to-apply 
nonlinear boundary conditions has been developed. This method may be useful 
in calculating flow about winged configurations. Finally, a simple integral 
boundary-layer method for compressible turbulent flow with a two-level eddy 
viscosity model has been developed. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 23, 1981 
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Experimental data (Klebanoff, ref. 13) 


^{l - erf[5(z/5 ~ 0.78)]} {Klebanoff, ref. 13) 
[l + 5.5(z/S) 6 ] (Cebeci and Smith, ref. 10) 


Present linear approximation 
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